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ABSTRACT 

IWGA Jl 346.5-6255 is a7?(9&4rX-ray source found within the radio lobes of the supernova remnant (SNR) 
G309.2-00.6. This source also appears to coincide with the bright and early-type star HD 119682, which is 
in the middle of the galactic open cluster NGC 5281. The radio morphology of the remnant, consisting of 
two brightened and distorted arcs of emission on opposite sides of the IWGA J1346. 5-6255 source and of a 
jet-like feature and break in the shell, led to the suggestion that IWGA J1346.5-6255/G309.2-00.6 is a young 
analog of the microquasar SS 433 powering the W50 nebula. This motivated us to study this source at X-ray 
and optical wavelengths. We here present new Chandra observations of IWGA J1346. 5-6255, archival XMM- 
Newton observations of G309.2-00.6, and optical spectroscopic observations of HD 1 19682, in order to search 
for X-ray jets from IWGA J1346. 5-6255, study its association with the SNR, and test for whether HD 1 19682 
represents its optical counterpart. We do not find evidence for jets from IWGA J1346. 5-6255 down to an 
unabsorbed flux of 2.6 x lO"'-' ergs cm"- s (0.5-7.5 keV), we rule out its association with G309.2-00.6, and we 
confirm that HD 119682 is its optical counterpart. We derive a distance of 1.2±0.3 kpc, which is consistent 
with the distance estimate to NGC 5281 (1.3 ±0.3 kpc), and much smaller than the distance derived to the SNR 
G309.2-00.6. We discuss the nature of the source, unveil that HD 1 19682 is a Be star and suggest it is a new 
member of the recently proposed group of 7-Cas analogs. The Chandra and XMM-Newton X-ray lightcurves 
show variability on timescales of hundreds of seconds, and the presence of a possible period of ^1500 s that 
could be the rotational period of an accreting neutron star or white dwarf in this 7-Cas analog. 

Subject headings: X-rays: individual: IWGA J1346.5-6255, IRXS J134633.6-625528, SNR G309.2-00.6 - 
stars: individual: HD 119682 - stars: emission-line. Be - open clusters and associations: 
individual: NGC 5281 - ISM: individual: G309.2-00.6 - ISM: abundances - 



1. INTRODUCTION 

X-ray observations of Supernova Remnants (SNRs) have 
been continuously shaping our understanding of the diversity 
of compact objects associated with them. While the Crab 
used to be the prototype example for the aftermath of a su- 
pernova explosion, there is now growing evidence for a di- 
versity of compact objects associated with SNRs. These in- 
clude the magnetically powered neutron stars or magnetars, 
the radio-quiet neutron stars, the high-magnetic field radio 
pulsars, and accretion-powered candidates. Out of the 236 
catalogued Galactic SNRs however, W50/SS 433 remains the 
only system confirmed to harbor a microquasar (a neutron star 
or black hole in a binary system emitting a two-sided semi- 
relativistic jet). The search for other such systems is needed 
not only to address the uniqueness of W50/SS 433, but also to 
study the connection between supernovae and microquasars. 
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To date, there are just six confirmed mi croquasars wi th mas- 
sive optical companions in our Galaxy ( Paredesll2005b . Since 
the physics of the accretion-ejection phenomenon in micro- 
quasars (and AGN) central engines is not well understood, 
each additional SNR-microquasar system helps greatly to 
constrain the theoretical models. Here we present a study 
of the G309.2-00.6/1WGA J1346.5-6255 system which has 
been proposed to be a young analog of the W50/SS 433 sys- 
tem. 

1.1. G309.2-00.6 and IWGA J1346.5-6255 

G309.2-00.6 was first identified as an SNR based on its 
nonthermal radio emission (Green 1974). A detailed study 
of the SNR was perfor med using the Austra lia Telescope 
Compact Array (ATCA, iGaensler et al] 119981) . The rem- 
nant has a distorted shell-like morphology (see Fig. [T] and 
Fig. |2]i with two brightened and distorted arcs of emission 
on opposite sides with diameters of 12' and 15', a jet-like 
feature and breaks in the shell. The SNR has a surface 
brightness of 5.4x10"^' W m'^ Hz"' Sr"' at 0.843 GHz 
dWhiteoak & Greenll 19961) . Interestingly, this morphology re- 
sembles almost exactly a scaled down version of the SNR 
W50 (see, e.g., Dubneret al. 1998), well known to be pow- 
ered by the microquasar SS 433. H I absorption measurements 
against the SNR yield a distance in the range 5.4-14.1 kpc, 
corre sponding to an age of (l-20)xl0^ yr (IGaensler et alj 
IT998I) . 

Infrared observations of the remnant show no o bvious emis- 
sion from G309.2-00.6 (IGaensler et all [T998h . but strong 
emission coincident with the nearby object RCW 80, iden- 
tified as an H II region based on its infrared and Ha emission 
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and flat radio spectrum. H I absorption measur ements indicate 
a dista nce of 5.4±1.6 kpc towards RCW 80. iGaensler et al.l 
(Il998l) suggest that the SNR may be interacting with the ad- 
jacent H II region, which would put the SNR at the lower 
range of the suggested distance and would imply an age 
<4,000 years. The striking resemblance of the radio morphol- 
ogy of G309.2-00.6 to W50 and the radio jet feature which 
aligns with the axis of symmetry defined by the two ears led 
IGaensler et al.l {l998) to suggest that the SNR is powered by 
a central source producing collimated outflows, or jets, in two 
opposed directions, like the microquasar SS 433. A central 
source, however, has not been detected in the radio band, with 
a 5ct upper limit of 0.4 mJy at 1.3 GHz. 

X-ray observations of G309.2-00.6 with the Advanced 
Satellite for Cosmology and Astrophysics (ASCA) have 
shown that the SNR is rich in therm al emission from a metal- 
rich ejecta-dominated plasma (Rak owski et al.ll2b01l) . indicat- 
ing a young SNR age (700^,000 years old). The X-ray spec- 
trum of the SNR is best fitted with a non-equilibrium ioniza- 
tion model with a column density A^h ^ 0.6 ± 0.3 x 10^^ cm"^, 
yielding an estimated distance of 4±2 kpc to the remnant. 

A bright ROSAT X-ray point-like source, 
IRXS J134633.6-6 25528 also known as 
IWGA J1346.5-6255 dWhite et alj [19961) . wifli a count 
rate of 0.06±0.02 counts s"' is located near the center of 
the radio emission from the SNR. This source is coincident 
with the early-type star HD 119682, the third brightest star 
in the middle of the galactic open cluster NGC 5281 (see 
^L2b . iRakowski et all (1200 Ih suggested the possibility that 
this source is related to the SNR based on its position and 
derived column density of (0.3-1.0)x 10^^ cm"^. 

1.2. HD 119682 

The star HD 119682 has been usually classified as one 
of the members of NGC 5281. Although it has been ob- 
served several times (see available V and B -V photome- 
try in Table [U no precise spectral type and luminosity class 
have ever been reported, preventing an accurate estimation 
of the distance to this star. HD 11968 2 was initially clas- 
sified as having a B spectral type by Cannon & Pic keririg 
([T9I8), and later on Stechenson & Sanduleak (1971) gave 
an OB+ spectral type. iMoffat & Vogti d 19731) suggest a lu- 
minosity class of III? and quote a distance of 1.30 kpc to 
NGC 5281 based on its unevolved stars. These authors also 
note that star #3, HD 119682, has an unusual position in 
the color-color diagram (CCD) and color-magnitude diagram 
(CMD ). Humphreys il975h reports a Be spectral type, while 
iFitzgerald et al.. (.197 9^) give a slightly different spectral type 
of 09e, although no luminosity class is specified. Based 
on the position of H D 119682 in the CC D and CMDs by 
iMoffat & Voi3 (Il973h . [Mermifliod (|1982h classified the ob- 
ject as a blue straggler. 

A detailed study of NGC 5281 was carried out by 
ISanner et al.l (1200 lb . They report further photometry for 
HD 1 19682 and explicitly state that they do not see evidence 
for a blue straggler nature in the CMDs when compared to 
other stars of NGC 528 1 . Based on isochrone fitting to the 
CMDs of the cluster members they provide the following pa- 
rameters: (m-M)o = 11.0 ±0.2 mag, d = 1580 ± 150 pc, 
Eb-v = 0.20 ± 0.02 mag, an age of f = 45 ± 10 Myr (or 
lo gf(yr) = 7.65 ±0.10) a nd a solar metaUicity of Z = 0.02. 

iLevenhagen & Leisted |2004) report the following physical 
parameters by fitting the spectral profile of HD 1 19682 with 
a non-LTE model: vsin/ = 220 ±20 km s"', Teff = 31910 ± 



550 K (or log reff(K) = 4.504±0.008), logg = 4.00±0. 10 dex, 
logf(yr) = 6.60±0.10 (orf = 4± 1 Myr), logL/Lg = 4.64 ± 
0. 10 and M = 18 ± 1 Mq. Note that the derived age is an or- 
der ofjMgnitude smaller than the one of NGC 528 1 reported 
by lSanner et alJ (1200 Ih . clearly reinforcing the blue straggler 
nature of the star. 

Finally, we note that near infrared photometry of 
HD 119682 is present in flie 2MASS catalog: J = 7.350 ± 
0.027, H = 7.184 ±0.033, and K, = 6.950 ± 0.023 mag 
(ICutri et alJl2003l:ISkrutskieetan i2006^. 

The goal of this paper is to study the nature of the 
IWGA J1346. 5-6255 source, its association with the optical 
star HD 1 19682, and the association of either source with the 
SNR G309.2-00.6. The previous proposal that G309.2-00.6 
represents a young analog of W50 motivated us to study 
the ROSAT X-ray source and any associated jets with the 
Chandra X-ray Observatory. We also complement our X-ray 
study with archival observations obtained with XMM-Newton. 
While a detailed study of the diffuse emission from the SNR 
is beyond the scope of the paper, we here focus on the possible 
association between the X-ray source and the SNR. We also 
present new optical observations of HD 1 1968 acquired with 
ESO's New Technology Telescope. These observations are 
targeted to pin down the spectral type and luminosity class of 
the star from which to obtain a good estimate of its distance. 

The paper is organized as follows. In ^ we summarize the 
observations undertaken at X-ray and optical wavelengths. In 
and|5] we present our Chandra, XMM-Newton, and the 
optical results, respectively. We finally discuss our results in 
^and summarize our conclusions in 

2. OBSERVATIONS AND DATA REDUCTION 

The ROSAT source IWGA J1346.5-6255 was observed 
with the Chandra X-Ray Observatory for 15.12 ks on 2004 
December 26, with the back-illuminated chip S3 of the Ad- 
vanced CCD Imaging Spectrometer (ACIS-S; G. Garmire*^), 
at a focal plane temperature of -120 °C ^. The data was pro- 
cessed using the standard CIAO 3.2.1 tools"^. Events with 
ASCA grades (0, 2, 3, 4, 6) were retained, and periods of 
high-background rates were removed. The net effective expo- 
sure time was 14.9 ks. The data was corrected for Charge 
Transfer Ineffic iency (CTI) using the Penn State corrector 
dTownslev et al.l 12002) and filtered using the filtering avail- 
able in acis_process_events (CIAO 3.2.1 with CALDB ver- 
sion 3.0.1). 

The IWGA J1346. 5-6255 source spectrum was created by 
selecting photons within a circle of radius 3 "5 and centered 
at the peak of the X-ray emission (see §3.1). The background 
was selected from an annulus centered on the X-ray source, 
with inner radius of 10" and outer radius of 15". The count 
rate for the source was 0.104 ±0.005 counts sec"' and the 
count rate for the background was 0.004 ± 0.002 counts sec"' . 
The spectrum was grouped to a minimum of 20 counts per 
bin and the energy range selected for spectral fitting was 
0.5-7.5 keV (the spectrum is background dominated above 
7.5 keV). The spectral analysis was performed using XSPEC 
version 1 1.2 

The SNR G309.2-00.6 was observed by the Euro- 
pean Photon Imaging Camera (EPIC) aboard XMM-Newton 



^ See http://cxc.harvard.edu/propos er/POG] 

' For a review of Chandra, see http://cxc.harvard.edu/cdo/about_chandra/overview-cxo.hti 

http://cxc.harvard.edu/ciao 
' ' http://xspec.gsfc.nasa.gov 
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dJansen et al.ll2001h for 40.6 ks on 2001 August 28 (obs ID 
0087940201). The EPIC-PN camera was operating in Ex- 
tended Full Frame mode, and the EPIC-MOS cameras in Full 
Frame mode. A thick filter was used to avoid problems with 
bright stars in the open cluster NGC 528 1 . The data were pro- 
cessed using version 6.0 of the XMM-Newton Science Anal- 
ysis Software (SAS). We created filtered event files using the 
EPlC-PN and MOS cameras in the 0.3-10 keV energy range. 
We then selected events for which the pattern was less than 
or equal to 4 (12) for PN (MOS), we removed the bad pixels, 
and cleaned the data by removing the times with high back- 
ground proton flares. The resulting effective exposure time is 
23.9 ks for PN, 27.8 ks for MOSl, and 27 ks for MOS2. The 
spectrum for IWGA J1346. 5-6255 was extracted using a cir- 
cle centered at the peak of X-ray emission: a = 13*" 46™ 32!^6, 
5 = -62° 55' 27" (J2000) and with a 0^5 radius. The spectrum 
of the diffuse emission from G309.2-00.6 was extracted us- 
ing a circle centered at the peak of the diffuse X-ray emission: 
a = O'^ 46™ 33^ S = -62° 50' 57" (J2000) and a radius of 
3fl25. This region is shown in Fig.|2] The background spec- 
trum for IWGA J1346. 5-6255 was extracted from a ring cen- 
tered at the source and extending from 0.'5 to l.'O; while the 
background spectrum for the diffuse emission was extracted 
from an annulus centered at the diffuse emission and of 5. '83 
and 9.'17 for the inner and outer radii, respectively. The ancil- 
lary and response matrix files were produced using the SAS 
commands arfgen and rmfgen, respectively. The spectra were 
subsequently grouped into a minimum of 50 counts per bin 
for the point source and 200 counts per bin for the diffuse 
emission. The spectral analysis was subsequently performed 
using XSPEC vll.2. 

We observed HD 119682 using the ESQ Multi-Mode In- 
strument (EMMl) on the 3.5-m New Technology Telescope 
(NTT) at La Silla, Chile, on 2003 June 5th. An Ha spec- 
trum was taken with the red arm in intermediate-resolution 
mode (REMD) and grating #6. The red arm is equipped with 
a mosaic of two thin, back-illuminated 2048 x 4096 MIT/LL 
CCDs and this configuration results in a nominal dispersion 
of 0.4 A pixel"' over the range AA 6440-7150 A. The reso- 
lution, measured on arc lines, is « 1.2 A. The blue spectrum 
was taken with the blue arm in intermediate-resolution mode 
(BLMD) and grating #12. The blue arm is equipped with 
a Texti-onik TK1034 thinned, back-illuminated 1024 x 1024 
CCD and this configuration results in a nominal dispersion of 
0.9 A pixel"' over the range AA 3820-4750 A. The resolution, 
measured on arc lines, is « 2.6 A. Image pre-processing was 
carried out with MIDAS software, while data r eduction was 
achieved with the Starlink packages CCDPACK dDraper et al.l 
I2OOQ) and FIGARO (Shorti-idge et al. 1993). 

3. C//AA'Di?A RESULTS 
3.1. Imaging 

In Fig.n we show the radio image of the SNR G309.2-00.6 
obtained with the Molonglo Observatory Synthesis Telescope 
(MOST) at a frequency of 0.843 GHz overlayed with the 
X-ray emission detected in the S3 chip of Chandra. The 
peak of the X-ray emission is located at a = O*" 46™ 32!'6, 
5 = -62° 55' 24" (J2000) with an error radius of 1"3 (90% 
confidence level). This source corresponds to the previously 
known /?(95Ar source IWGA J1346.5-6255. 

Our primary goal is to search for any jets associated with 
this X-ray source. Since jets from compact objects gener- 
ally peak at different X-ray energies than their powering en- 



gines, we created energy images for IWGA J1346. 5-6255 
in the soft (0.3-1.3 keV), medium (1.3-2.2 keV), and hai'd 
(2.2-10 keV) band. The energy boundaries were chosen as 
to obtain a similar number of counts (^-^500) in each image. 
The images were then smoothed using the CIAO tool csmooth 
with a circular Gaussian kernel and are shown in Fig. [3] The 
initial smoothing scale was 1 pixel. The smoothing scale 
was increased (to a maximum 2 pixels) until the number of 
counts under the kernel exceeded a signal-to-noise ratio of 2. 
Where the signal-to-noise ratio was less than 2, the image was 
smoothed on a scale of 2 pixels. The axis overlayed on top of 
the images, and running from northeast to southwest, delin- 
eates the axis of symmetry of the radio lobes of G309. 2-00.6 , 
~ 45° counter-clockwise from north jGaensler et al.l [T998I) . 
Interestingly, there is a hint of extension in the soft energy 
band along this axis. However, any extension observed in 
Fig. |3] (purple, dark blue) has a count rate of less than 0.7 
counts pixel"', indicating that the apparent extension is most 
likely background noise. 

In order to further determine whether IWGA J 1346.5-6255 
is point-like or extended, we subsequently compared its 
spatial characteristics with Chandra's, point spread function 
(PSF). We generated a PSF image at an offset axis of 0.'48 
(the location of the source on the S3 chip) and at an energy 
of 1.5 keV (characteristic of the source's energy histogram). 
We subsequently normalized the PSF image to the source 
counts (^-^1450 counts in the 0.3-10 keV range) and used 
the PSF as a convolution kernel when fitting the source. Fit- 
ting a two-dimensional Gaussian distribution to the data, we 
found that the full width half maximum (FWHM) was only 
slightly larger than the FWHM of the PSF which is assumed 
to be circular (1"02 and 0."73, respectively). The data and 
PSF are shown in Fig. |4] Since the difference is less than 
one pixel (0"4920 ± 0"0001), the data is consistent with a 
point source. A similar profile was created for the soft (0.3- 
2.4 keV) and hard energy band (2.4-10 keV), by evaluating 
the PSF at 1.1 keV and 2.5 keV, respectively. The FWHM of 
the data and PSF in the soft band was found to be 0."97 and 
0"73, respectively. For the hard energy band, the FWHM of 
the data and PSF was 1"09 and 0"77, respectively. We con- 
clude that the data in either band are consistent with a point 
source, and therefore rule out the presence of jets associated 
with IWGA J1346. 5-6255, down to an unabsorbed flux of 
2.6 xlO"'"* ergs cm"^ s"' in the 0.5-7.5 keV range (see next 
section). 

3.2. Spectroscopy 

In order to address the nature of IWGA J1346. 5-6255, 
we fitted its spectrum with thermal and non-thermal mod- 
els, modified by interstellar absorption. In Table |2l we sum- 
marize the parameters of the best fit models. Errors are 
at the 90% confidence level throughout the paper, unless 
otherwise mentioned. We find that blackbody and thermal 
bremsstrahlung models do not provide acceptable fits. The 
absorbed power-law model provides an adequate fit with the 
following parameters: a column density A^h = 0.18;!;Qjj7 x 
10^~ cm"^, a photon index F = 1.4 ±0.2 and an unabsorbed 
flux of l.lil^ jX 10"'^ ergs cm"^ s"' in the 0.5-7.5 keV range 
(reduced x^=0.46, i^=61 degrees of freedom). The spectrum 
and fitted model are shown in Fig. |5] while in Fig. |6] we 
show the 68, 95 and 99.7% confidence contours for A^h and 
r. By examining the residuals in Fig. |5] some line emission 
appears to exist around 1 keV. For this purpose, we added 
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a MEKAL component (which would, e.g., account for any 
emis sion associated with coronal emission from a normal star, 
see ^6. 51 for more details). We find that the fit improves sig- 
nificantly (xl, = 0.341, h' = 59), with an F-test probability 
of 4.5x10"^ that this improvement occurs by chance. This 
best fit model is shown in Fig. [T] and the parameters are: A^h 
= 0.19!°;,')^ X 10^2 cm-2, T = 1.24+0;2^, kT = 0.97!°;^ keV, 
with an unabsorbed flux in the 0.5-7.5 keV energy range of 
l-18;!;o29^ 10 '^ ^^S,^ cm"^ s"', very similar to the one ob- 
tained with a pure absorbed power law. 

In order to test whether the X-ray emission is coronal emis- 
sion from a normal star (see %.5\ . pure MEKAL models 
were attempted. A single temperature MEKAL model yields 
a statistically acceptable fit, however the temperature is high 
and unconstrained (see Table |2]i. A two-component MEKAL 
model yields as good a fit as the power-lawH-MEKAL model 
but with the temperature of the hard component poorly con- 
strained: A^H = X 10^2 cm-2, kTi = 0.95!{J:|^ keV; 
kT2 = 80 (> 14) keV, and a reduced xl = 0.342 for = 59 de- 
grees of freedom. The temperature of the hard component, of 
~ 9 X 10*^ K, is unrealistically high and poorly constrained- 
we show in the next section that the XMM-Newton data helps 
better constrain this component. We conclude that while the 
two-component MEKAL model and the power-lawH-MEKAL 
model both yield statistically acceptable fits, the Chandra 
spectrum is best fitted and its parameters better constrained 
using the power-lawH-MEKAL model. 

In order to quantify the upper limit on the flux from any 
unseen jets associated with IWGA J1346. 5-6255, we added 
a power-law component to the best fit power-lawH-MEKAL 
model above. We then froze the parameters of the model 
for IWGA J 1346.5-6255 and allow the photon index for 
any 'unseen' jets to vary between 1.6-2.0, which is reason- 
able based on observations of other jet sources. We find 
that the corresponding upper limit on the unabsorbed flux is 
2.6 X 10"'"^ ergs cm~^ s"' in the 0.5-7.5 keV range. At a dis- 
tance of 1.3 kpc (see §6), this translates to a luminosity of 
~5.3x 10-^' ergs s"'. 

3.3. Timing 

To check for variability of IWGA J1346. 5-6255, we cre- 
ated lightcurves with Chandra data binned at 100, 120 and 
150 s time intervals. The obtained lightcurves display vari- 
ability with an amplitude at the ±30% level {la value) on a 
timescale of a few 100 seconds. 

We have searched for periodic signals using standard 
techniques like th e Phase Dispersion Minimi zation (PPM, 
StellingwerHl 19781) and the CLEAN algorithm ([Roberts et al.1 
19871) . The PDM periodograms obtained for the Chandra 
data, binned at either 100, 120 or 150 s, reveal a periodicity 
of 1500!5[f s. The obtained PDM statistic 6 is in the range 
0.89-0.92. An F-test reveals that this periodicity could be the 
result of a random fluctuati on with a probability of ~40% (see 
details in lStellingwerlll97 8). This is not strange given the fact 
that the noise witihin individual PDM phase bins is two times 
higher than the overall detected variability. We note that sim- 
ulations with even lower signal-to-noise ratios have revealed 
the good perf ormance of PDM (as can b e seen for the primary 
periodicity in 'Otazu et 111120021 12004|) . On the other hand, 
two additional possible periodicities of ^3050 and ^4600 s 
are detected with lower significance. Given their spectral pro- 
files, it is clear that they are subh armonics of the ^^1500 s 
periodicity (see IStelUngwergi 19781 for details). The CLEAN 



algorithm displays less significant maxima at ^^1550 s, possi- 
bly because the signal is not perfectly sinusoidal but is better 
accounted for by a saw-tooth pattern with a slow increase and 
a relatively fast decay. However, sinusoidal fits to the binned 
data provide a possible period in the range 1545-1555 s, with 
a variability amplitude of 16%. We show in Fig.[T4]the back- 
ground subtracted Chandra lightcurve binned at 150 s inter- 
vals and normalized to its average value. We also present 
the same data smoothed with a running window of 450 s to 
enhance the 1500 s modulation, and a sinusoidal fit to the 
binned data, which approximately follows the relative max- 
ima and minima of the smoothed data. In §6.5, we discuss 
these results, together with the light curves obtained using the 
XMM-Newton data (next section). 

4. XMM-WEWTOA' RESULTS 

4.L IWGA J 1346.5-625 5 

In order to fit the XMM-Newton spectrum of 
IWGA J1346.5-6255, we attempted single component 
(thermal and non-thermal) and two-component models, 
as we did above for the Chandra data 03.2l i. The XMM 
data allowed us however to extend the fits up to 8.5 keV. 
As with Chandra, we found that the power-law model 
provides an adequate fit with A^h = 0.22 ±0.02 x 10^- cm"^, 
a photon index r=1.69;!^ jjg, and an unabsorbed X-ray flux 
of 2.22 ± 0.02 X 10-'^ ergs cm'^ s"' in the 0.5-8.5 keV 
energy range, with a reduced x?,= l-04 (for t^=418 degrees of 
freedom). As with the Chandra section, errors are at the 90% 
confidence level, unless otherwise mentioned. The spectrum 
is shown in Fig. |8] and the confidence contours for A^h and F 
are shown in Fig.|9] 

By comparing the XMM-Newton to the Chandra param- 
eters, we note that the XMM-derived A^h and F values are 
slightly higher than those derived from the Chandra analy- 
sis 03.2| i. The contour plots however (Figs. |6]and|9]l show 
that within errors, the values are consistent with each other, 
with the XMM-Newton values being much better constrained. 

We subsequently tested for whether pile-up affected the 
Chandra spectral analysis (the XMM-Newton spectrum is not 
affected by pile-up). Using Webpimms and a count rate of 
^0.10 counts s"' in the 0.5-7.5 keV, we estimate a pile-up 
fraction of <12% in ACIS-S3. We then corrected for pile- 
up using the pileup model in XSPEC. We find that the A^h 
and F values increase slightly and are closer to the best val- 
ues obtained with XMM-Newton, but the range is still con- 
sistent with our previous Chandra estimates and the contour 
levels are very similar to those shown in Fig. |6] After cor- 
rection for pile-up, we find that A^h = 0.20:I:o j], x 10^^ cm"^, 
and F = 1.53;!;o33. corresponding unabsorbed flux is 

2 ^ 10"'^ ergs cm"^ s"\ consistent with the value be- 
fore pile-up correction, but only accounting for ~50-55% of 
the flux detected by XMM-Newton in the 0.5-7.5 keV range. 
We conclude that while the spectral parameters are not sig- 
nificantly affected by pile-up, the difference in the flux esti- 
mate between the Chandra and XMM observations may be 
partly attributed to some contamination by the background 
and likely evidence for variability in the source (see ^ for 
a discussion on the variability of the source). 

We then attempted pure MEKAL models. A one- 
temperature component MEKAL model provides a statisti- 
cally acceptable fit, however the model overestimates the 
X-ray emission between 6-7 keV and is therefore dis- 
carded. Just like with the Chandra data, a two-temperature 
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MEKAL model provides a better fit with the following pa- 
rameters: A^H= 0.16±0.01 xl022 kTi= 1.7±0.3 keV, 
kT2= 13.0!|^ keV, xl = 1-025 (i/=416 degrees of freedom). 
In Fig. [To] we show the EPIC spectra fitted with this two- 
component model and in Table[3] we summarize the results in- 
dicating that the harder component dominates the X-ray emis- 
sion. In Fig. [n] we show the 68, 95 and 99.7% confidence 
contours for the fitted temperatures of the two-component 
MEKAL model. Note that the harder component's temper- 
ature is much better constrained with XMM-Newton than with 
Chandra. This temperature is still high for coronal emission 
but not inconsistent with the ones measured in the new 7 Cas- 
like objects (see the discussion section). 

Lastly, a two-component power-lawH-MEKAL model pro- 
vides the best fit with the lowest x^, with the following pa- 
rameters: A^H = 0.2l!°;U^xl022 cm-2, kT = 1.4!|jj keV, T = 
1.60±0.07, reduced x^=1.005 (z/=416 degrees of freedom). 
The corresponding fit is shown in Fig.[T2]and the parameters 
summarized in Table [3] As in the MEKALh-MEKAL fit, the 
harder component (now the power-law) dominates the X-ray 
emission. 

Finally, it is worth noting that the power-law fit (see Fig.|8]l 
indicates that a bump near 6-7 keV is unaccounted for. 
Adding a Gaussian line to the power-law component im- 
proves the fit (x4j5= 1.012), with an F-test probability of 0.1 17 
that this improvement occurs by chance. The line energy 
parameters are: E\i„s = 6.77;I^ [2 keV, criine=0.19^{J [g4 keV, 
and an unabsorbed line flux of 4.8 x 10"'"* ergs cm"^ s"'. 
Adding a Gaussian line near 6.7 keV to the two-component 
(MEKALh-MEKAL or MEKAL-npower-law) models does not 
improve the fits significantly. 

We conclude that the two-component MEKALh-MEKAL, 
power-law+MEKAL, and power-law+Gaussian line, all pro- 
vide good fits to the XMM-Newton data, although the power- 
law+MEKAL provides the best fit, a result that is consistent 
with the Chandra study. 

As for the Chandra data, we also obtained background sub- 
tracted lightcurves binned at 150 s and normalized to their av- 
erage values. As can be seen in Fig. [15] IWGA J1346. 5-6255 
was significantly variable, with an amplitude at the ±45-55% 
level (1(7 values) on timescales of a few 100 seconds. A 
timing analysis reveals different possible signals at ^-^4570, 
^8300 and ^12000 s in the different lightcurves and us- 
ing PDM or CLEAN. However, sinusoidal fits to the data 
with such periods do not converge. Since a marginal sig- 
nal is nearly always detected at ^^1500 s with PDM and/or 
CLEAN, we have tried a sinusoidal fit around this periodic- 
ity, and found that convergence is achieved for a period of 
1485±10 s and ampHtudes of 10, 20 and 16% in the PN, 
MOSl and MOS2 lightcurves (to be compared to the 16% 
amplitude in the Chandra data, §3.3). We have also analysed 
a combined data set of all normalized XMM data, and found 
similar results: the three possible periods reported above and 
a more clear signal than in individual data sets at 1486±2 s 
(and 2 subharmonics for PDM), with an amplitude of vari- 
ability at the 13% level, hence closely matching the Chandra 
results. An F-test to the PDM results reveals a probability 
of up to 60% of this signal as being due to a random fluc- 
tuation. Although from the statistical point of view we can 
not argue that the periodicity exists, the detection of 2 subhar- 
monics and the fact that we detect a very similar periodicity 
with Chandra and XMM data gives us confidence that the pe- 
riod, or quasi-period, is real. Future observations can probe 



this issue. 

4.2. G309.2-00.6 

We subsequently fitted the spectrum of the diffuse emission 
from the SNR G309.2-00.6. While a detailed study of the 
SNR is beyond the scope of this paper, we here present the 
spectral fit in order to compare its column density to that of 
IWGA J1346.5-6255 (see Fig. |2]i, and therefore test for the 
association between the compact source and the SNR. Since 
G309.2-00.6 is beUeved to be a young SNR (Rak owski eTal] 
120011) . we fitted its spectrum with a non-equilibrium ioniza- 
tion (NEl) model with variable abundances (vnei in XSPEC, 
lBorkowskiet"ani2001h . A vnei model with solar abundances 
did not yield an acceptable fit (xi = 10.5, i' = 169 degrees of 
freedom ), a result that is consi stent with the previous ASCA 
study bv lRakowski et al.] (120011 ). We subsequently varied the 
abundances of O, Ne, Mg, Si, S, Ar, Ca, Fe and Ni; and found 
that the model requires an over-solar abundance of Si and S, 
and an under-solar abundance of O, Ne, Mg, Ca and Fe. The 
plasma temperature is kT = 2.0±0.6 keV and the column den- 
sity is A^H = 0.65;!;o25 xlO^^ cm"^ (errors here are at the 3a 
level, see Fig. \T3[ . This best fit A^h value is at least three 
times higher than that derived for IWGA J1346. 5-6255 (see 
Tables [2] and [3). Even considering the upper 3 a value for 
IWGA J1346.5-6255 with Chandra or XMM and the lower 
3(7 value for G309.2-00.6, A^h of the diffuse emission from 
the SNR is still higher. Therefore, based on the A^h values 
alone, we can conclude that IWGA J1346. 5-6255 is unlikely 
associated with the SNR. We note that this result does not 
agree with Rakowski et al.'s result, likely due to the lack of 
resolving power with ASCA which did not allow a study of the 
SNR independently of the IWGA J 1 346.5-6255 point source. 

5. OPTICAL RESULTS 

The normalized classification spectrum of HD 119682 is 
shown in Fig. [16] It is covered in emission lines, most cor- 
responding to Fe II transitions, and looks typical of Be stars. 
Spectral classification is complicated by the presence of many 
weak emission lines which render identification of weak ab- 
sorption features very difficult. The spectral type can be es- 
timated from the presence of both Si III and Si IV lines and a 
very weak He II 4686 A line to be approximately BO. 5. The 
determination of the luminosity class is complicated by the 
emission components and broadened lines, but the weakness 
of Si IV 4089 A indicates that the star is not very luminous and 
suggests a luminosity class V. This is reinforced by the value 
of log ^ = 4.00 ± 0.10 dex found by 'L evenhagen & Leisted 
(120041) . which is consistent with the log § = 3.9 values for main 
sequence stars but clearly incompatible with the \ogg =3.5 
values for giants listed in Martins et al. ( 2005) . Therefore, we 
will adopt a spectral classification of BO. 5 Ve. 

We note that the C III 4650 A line is abnormally weak 
for a star of this spectral type, as it is always stronger than 
O II 4642 A for any star eai-lier than B2. The N II 3995, 
4041 A lines are moderately stronger than expected, suggest- 
ing that there is a certain degree of CNO anomaly. This could 
explain the abnormal C III 4650 A/O II 4642 A ratio, as C 
could be moderately depleted and the A 4642 A line could 
have a contribution from N II 4640 A. 

The Equivalent Width of Ha emission line is -24 ± 2 A. 
The shape is typical of a Be star with strong emission, and 
shows no discernible structure or asymmetry. 
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6. DISCUSSION 

6.1. Astmmetry ofHD 119682 and IWGA J1346.5-6255 

The most precise position and pro per motion of HD 1 1968 
are given in IZacharias et al.l (12004 ^. Considering that our 
Chandra observations were conducted on year 2004.987, 
we obtain the following J2000 coordinates for that epoch: 
a = B** 46™ 32?579 ± 0?005, S = -62° 55' 24"17 ± 0."03. 
This is perfectly consistent with the Chandra position of 
IWGA J1346.5-6255: a = 13^^ 46" 32?6, 5 = -62° 55' 24" 
(with 1"3 error radius at 90% confidence level). However, 
given the high density of stars in the region where the X-ray 
source is located, we have estimated the chance coincidence 
probability of a 2MASS object (assuming negligible error in 
position, typically < 0."1) falling within the Chandra error 
circle in position. We have done this by counting the 2MASS 
objects located within circles of different radius (1, 2, 3, 5, and 
10') around the Chandra position, and found that the proba- 
bility is always smaller than 4% for objects brighter than those 
present in the 2MASS catalog for this field: 7 ~ 18, // ~ 17, 
and Ks ~ 16 mag. In fact, the absolute difference in position 
is just 0"22, much smaller than the 1 "3 Chandra position er- 
ror, providing a chance coincidence probability of ~ 0.1%, 
which practically confirms without any doubt the association 
between HD 11968 and IWGA J1346.5-6255. 

6.2. Distances to HD 119682/lWGA J1346.5-6255 and SNR 
G309.2-00.6 

Our spectral type and luminosity class estimates of 
HD 1 19682 allow us to compute the distance to the star, and 
hence to IWGA J1346.5-6255. For a B0.5 V star the ab- 
solute mag nitude is My = -3.6 ± 0.4 mag (slightly extrapo- 
lating from iMartins et al.l |2005|) while the intrinsic color is 
(B-V\) = -0.28 ± 0.02 dJolmsonllilli). The average of the 
photometry in the literature, reported in ^ is V = 7 .97 ± 0.06 
and B-V = 0.16 ±0.05 mag. The color excess is then Eb-v = 
0.44 ± 0.07 mag and the visual extinction is Ay = 1.41 ± 
0.24 mag(withAy = (3.30 +0.28(fi-y)o+0.04£jj-v)£(g-V) 
from lSchmidt-KaleJ 1 9825 . We finally derive a distance mod- 
ulus of V -My -Ay = 10.2 ± 0.5 and a distance to the source 
of t/= 1080±230pc. 

However, HD 1 19682 is in fact a Be star, so the contribution 
from the stell ar envelope to V and Eg-y should be considered. 
According to iDachs et"aLl ( Il988l) and using our measure of 
the Ha Equivalent Width (EW) of -24 ±2 A, we should sub- 
tract a contribution from the envelope of +0.05 ± 0.02 mag 
to Es-y and -0.08 ± 0.02 mag to V, leading to star values 
of y = 8.05 ± 0.06 and Eg-y = 0.39 ± 0.07 mag. This pro- 
vides a visual absorption of Ay = 1.26 ±0.23 mag, a dis- 
tance modulus of V -My-Ay = 10.4 ±0.5, and a distance 
of d = 1200 ± 260 pc. We note that our measure of the EW of 
Ha is not simultaneous to any of the used photometric data. 
However, since we use the average photometry, we expect 
these corrections to be approximately valid even if using a 
single value of EW. This distance estimate to HD 1 19682, of 
1 .2 ± 0.3 kpc, is in agreement, within errors, with the distance 
estimates to NGC 528 1 reported in the literature, the average 
of which is 1.3 ±0.3 kpc. Therefore, this will be the distance 
to HD 119682 adopted hereafter. We note that this is much 
lower than the distance estimates to SNR G309. 2-00.6 ap- 
pearing in the li terature: 5 .4 ± 1 .6 k pc (Gae nsler et al.lfT998h 
and 4 ± 2 kpc ( Rakowski e"taDl2001h . 

Finally, by using Ay = 1.26 ± 0.23 mag and the relationship 
of iPredehl & Schmitti ([T995h we obtain a hydrogen column 



density of A^h = (2.3 ± 0.4) x 10^' atoms cm"-. This is per- 
fectly consistent, within errors, with the values obtained from 
the X-ray spectral fitting of the Chandra and XMM-Newton 
data (see Figs. |6] and |9] and Tables|2]and[3]l. Therefore, there 
is no need for intrinsic absorption around the X-ray emitting 
zone. 

We can now estimate the distances to the X-ray sources 
based on the measurements of the hydrogen col umn density, 
their c onversion to extinction according to Predeh l & Schmitg 
( 119951) ( < Nn/EB-y > = 5.55 x lO^' cm'^ mag"'), and the 
value of the extinction per unit distance in the d irection of th e 
sources, which is of EB-y/d ~ 0.3 mag kpc"' (lLuckelll978h . 
For this purpose we use the most accurate values of A^h, ob- 
tained with the XMM-Newton data. For 1 WGA J1346.5-6255 
we consider the power-law model fit with A^h (2.2 ± 0.4) x 
10^' atoms cm"^ (3a uncertainties) and obtain a distance 
of c/ = 1.3 ±0.3 kpc, fully consistent with the distance to 
HD 119682 and NGC 5281. (We note that the column den- 
sity derived for the best fit power-law+MEKAL or power- 
law+Gaussian is close enough to that derived using the power- 
law and thus yields a similar distance estimate within error). 
For the diffuse emission of SNR G309.2-00.6 we use the vnei 
model fit with A^h = 6.5;!;2 5 x 10^' atoms cm"^ (3a uncertain- 
ties) and obtain a distance of d < 3.9t^i kpc, consiste nt with 
the 4 ± 2 kpc value obtained by lRakowski et al.] (12001) from 
their ASCA fit to the SNR. This distance is clearly incompat- 
ible with the distance to HD 119682/lWGA J1346.5-6255, 
but consistent within errors with the dis tance of 5 .4 ± 1 .6 kpc 
to the RCW 80 H II region proposed by lGaensler et al.1 (Il998h 
to be associated with SNR G309.2-00.6. 

6.3. Circumstellar envelope 

The intrins ic (J-Ks)o infrared co lor of a normal B0,1 V 
star is -0.22 (iRuelas-Maivorgal 1 99 lb . The 2MASS photome- 
try provides a value of J-Ks = 0.40 ± 0.04, leading to a total 
infrared color excess of Ej-k^ = 0.62 ±0.04. On the other 
hand, from the visual absorption of interstellar origin found 
above. Ay, is = 1 .26 ±0.23 mag, and using the relationships by 
iRieke & Lebofekvl (fT985h we found interstellar infrared ab- 
sorptions of Ay. IS = 0.36 ±0.06 and Ajf ,, is = 0.14 ±0.03, and 
finally Ej-k^^ is = 0.22 ±0.07. Therefore, in addition to the in- 
terstellar color excess, there is an extra contribution to Ej-k, 
of 0.40±0.08. This is typical of Be stars with strong emission 
lines, usually adscribed to their circumstellar envelope. 

6.4. Age ofHD 119682 

iLevenhagen & Leistej (l2004h report an age of 4 ± 1 Myr 
and a mass of M = 18 ± 1 M^n for HP 1196 82 based on 
the evolutionary tracks by ISchaller et al.l (Il992h. In fact, ac- 
cording to the new evolutionary tracks bv iMevnet & Maeded 
(2003) (see their Fig. 5), the luminosity a nd eff ective temper- 
ature reported by ILevenhagen & Leisterl (2004) are in agree- 
ment with a 1 8 Mq star at the beginning of the main sequence 
or with a slightly evolved 20 Mp, star. However , since new 
calibrations such as the one of' Martins et al.ll2005l favor lower 
masses, we will consider as a strict upper limits for the age 
of HD 119682 the main sequence lifetime of a 15 Mq star 
On the other hand, since HD 119682 has a high rotational 
velocity of vsin; = 220 ±20 km s~ ' we w ill use the evolu- 
tionary tracks bv lMevnet & Maeded (l2003h for an initial rota- 
tional velocity of 300 km s~\ leading to a total main sequence 
lifetime of 14.5 Myr (see their Table 1). Therefore, even if 
assuming the lowest possible mass and a high rotational ve- 
locity, the lifetime in the main sequence is less than 15 Myr, 
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i.e., 3 times lower than the estimated age of NGC 528 1 by 
ISanner et al.l (1200 Ih . This clearly confirms the blue straggler 
nature o f HP 1 19682 , which is not clear in the V vs. {B-V) 
CMD of ISanner et all (1200 Ih . slightly affected by the redden- 
ing of the envelope, but cl early evident wh en looking at the V 
vs. (U-B) CMD of Moffat & Vogj (UqtI . 

6.5. On the nature ofHD 119682/lWGA J1346.5-6255 

The unabsorbed X-ray luminosity of IWGA J1346.5-6255 
from our Chandra observations is L(o.s-7.5 keV) = 
2.5!o:^(c//1.3 kpc)2 X 10^^ ergs s"' (using the range of 
fluxes tabulated in Table |2]l. From the XMM-Newton 
observations and using the range of fluxes obtained for 
the various models (Table |3j, we obtain a luminosity of 
^^(0.5-8.5 keV) = ^■5^ij{d/\3 kpc)^ X lO"*^ crgs s"', which cor- 
responds to L(o.5-7.5 keV) = ^ ■'^'^liA / 1 3 kpc)^ X 10"^^ ergs s"', 
about a factor of 1.5-2 higher than the Chandra value, 
regardless of the model used. In any case, this is clearly 
lower than that of Be/X-ray binaries containing neutron stars, 
eve n the faint persi stent ones (which have Lx ~ lO''^ ergs s~', 
see iReig & Rochd |i999). Moreov er, no X-ray pu l sations 
have ever been detected, although iRakowski et al.l (1200 Ih 
could only constrain the pulsed fraction to be less than 
-85% in the frequency range 0.002-32 Hz (0.03-500 s) at a 
99.99% confidence level. Therefore, the currently available 
information does not favor the hypothesis of a binary system 
with an accreting neutron star as the compact object. 

Another possibility would be to have a relatively quies- 
cent accreting black hole in a binary system. Indeed, the 
photon index of IWGA J1346. 5-6255 is similar to those 
foun d in accreting blac k holes while in the lo w/hard state 

i see [ Fender et al.l 1200 41 and references therein). iGallo et akl 
2001 " found a correlation between the 2-11 keV X-ray 
flux and the cm radio flux density for black holes in the 
low/hard state. >From our power-law model fit to the XMM- 
Newton data, we obtain an unabsorbed flux of F2-U keV = 
1.7 X 10"'^ ergs cm"^ s"' and compute an expected radio flux 
density of 0.29 ± 0.20 mJy at a distance of 1 .3 kpc. Given the 
uncertainties, such a source could be above or below the 3cr 
level in the ATCA image of Gaensler et al. ( 1998). Deeper ra- 
dio observations, with a noise smaller than 0.02 mJy, and with 
long baselines to avoid contamination from the SNR, would 
be needed to properly search for a radio counterpart. Only 
simultaneous X-ray observations could be useful to exclude 
the black hole scenario, although this would also be uncertain 
due to the scatter of the radio/X-ray correlation. However, 
we note that there is no clear evidence that a Be-nblack hole 
binary system has been found up to now. In addition, the X- 
ray spectra are better fitted including one or several MEKAL 
models, which is not the case for black holes in the low/hard 
state. Therefore, the scenario of an accreting black hole orbit- 
ing the Be star is not favored by the available data. 

OB stars are known to be X-ray sourc es, presumab ly due to 
shocks arising in their stellar winds (see Giidel '2004' for a re- 
view). Based on Einstein data, Pallavicini et al. ( 1981) found 
that the X-ray luminosities of —30 OB stars followed approx- 
imately the empirical law Lx — (1 .4±0.3) x 10~^Lboi- A more 
complete stu dy based on ROSAT da ta for more than 200 OB 
stars allowed feerghofer et al.l d 19971) to propose the existence 
of a correlation with two different power laws above and be- 
low Lboi = 10^^ ergs s"'. The ^ o fX/L p, = 4.64 ± 0.10 value 
given by iLevenhagen & Leister! (|2004|) provides a star bolo- 
metric luminosity of Lboi = ( 1 .7 ± 0.4) x 10^^ ergs s~' , leading. 



according to iBerghofer et al.l (1 1997b . to an X-ray luminosity 
of Lx = (2;!;j) X 10^' ergs s~'. In contrast, the unabsorbed 
power-law X-ray luminosity from our Chandra observations, 
i(0.5-7.5 keV) = 2.5^q jid /13 kpc)^ X lO-*^ ergs s"' is around 
one order of magnitude higher, while the XMM-Newton lu- 
minosity is even a factor of about 2 higher. In fact, the X-ray 
emission of HD 1 1 9682 is not only stronger than expected, but 
obviously harder than that of isolated stars: our XMM-Newton 
two-component MEKAL fit to the data requires a high tem- 
perature of kT = 13.0!|5 keV (14.3!|^ x 10^ K), more than 
an order of magnitude higher than usual in main sequence B 
stars, and a factor of at le ast four higher than the most extreme 
cases (ICohen etalJll997h . 

As a matter of fact, the luminosity and spectral shape of 
HD 119682/lWGA J1346.5-6255 are similar to those dis- 
played by the notoriously peculiar Be star 7 Cas: L(2-io keV) = 
6 X 10^^ ergs s"' and two-component MEKAL spectrum 
wifli temperatures of 0.05 ± .01 keV and 12 3 ± 0.6 keV 
(IOwensetalJll999h . Recendy, ISmith & Baional (|2006) have 
proposed that the star HD 1 10432, which shows L(2-io keV) — 
5 X lO-*^ ergs s~' and a MEKAL spectrum with kT = 10.6 ± 
1.9 kev (lTorrei6n&Onil2001h . is an almost perfect twin 
of 7 Cas and have mentioned the possibility of the exis- 
tence of a class of objects displaying similar ch aracteristics, 
which they dub "7 Cas analogs". In addition, iMotch et alJ 
(l2005h have very recently presented a summary of X-ray 
and optical properties of their proposed currently existing 
five "7 Cas-like objects" (plus 7 Cas itself). The properties 
of HD 119682/lWGA J1346.5-6255, to be compared with 
their Table 1, are: spectral type B0.5 Ve, Ha EW of -24 A, 
Lx (0.2-12 keV) = 3.6-6.6 X 10^^^ ergs s~' (from the Chandra and 
XMM-Newton power-law model fits, respectively), kTsoh = 
1.7 keV and kT^^i-i = 13 keV, and F —1.7. These properties 
are practically identical to those of SS 397, with a slightly 
different equivalen t width of Ha. Furth ermore, one of the ob- 
jects in the list of iMotch et al.l ( l2005h is a blue straggler in 
the 50 Myr old open cluster NGC 6649, closely matching the 
observed properties of HD 1 19682, also a blue straggler and 
located in the -45 Myr old open cluster NGC 5281. 

In the previous spectroscopic study, we showed that the X- 
ray flux of the IWGA J1346.5-6255 source is a factor of 
—2 brighter with the XMM observation than with the Chan- 
dra one, obtained 3 years later. This could be attributed 
to intrinsic variability of the source. The short-term vari- 
ability detected on timescales of a few hundred seconds is 
simila r to that seen in other 7 Cas-like objects CMotch et alj 
l2005h . Moreover, the periodicity of —1500 s we have de- 
tected in HD 119682/lWGA J1346.5-6255, with a pulse 
fraction of 13-16%, is reminiscent to, although not so clear 
than, the oscillation found in the 7 Cas analog HD 161103, 
with a period of 3250 ± 350 s and a pulse fraction of 24% 
(iLopes de Oliveira et al.l |B006). However, in the case of 
HD 119682/lWGA J1346.5-6255, this possible period has 
been detected in two different observations performed with 
two different satellites, indicating that this signal is probably 
stable along time (although other superimposed signals are 
detected in the XMM-Newton data set). 

In summary, HD 119682/lWGA J1346.5-6255 shares the 
following properties with the other 7 Cas analogs: spec- 
tral fits including thermal components are better than simple 
power-law fits (see Tables |2] and [3]); in the two-temperature 
MEKAL model fits, the cooler component is much fainter 
than the hard one (see Tables |2] and [3]); the total photoelec- 
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trie absorption in X-rays is not very different from that due to 
the interstellar medium; the X-ray luminosity is in the range 
Lx (0.2-12 keV) = lO-'^-lQ-'-' ergs s"'; there are no reported X- 
ray outbursts; the flux varies by a factor of a few from one 
observation to the other; variability on short timescales is de- 
tected, showing different quasi-periodic signals; the presence 
of an iron line near ~6.7 keV (see §4.1 and Fig.|8]l, difficult 
to interpret in the black hole scenario; a moderate N enhance- 
ment; and a circumstellar envelope that contributes to the NIR 
reddening. 

Therefore, we conclude that the overall properties of 
HD 119682 indicate that it is most likely a new 7 Cas ana- 
log. 

The nature of 7 Ca s analogs is still a matter of debate (see 
iMotch et al.ll2005l and rLopes de Oliveira et al. ''2006 for recent 
discussions). Possible scenarios to explain the properties of 
the detected X-ray emission are the magnetic interaction be- 
tween the star and its circumstellar decretion disk, or accre- 
tion on to a compact object (neutron star, white dwarf, or 
even a black hole). We note that the possible ^^1500 s period 
we have detected in HD 1 19682/1 WGA J1346.5-6255 could 
be the rotational period in the neutron star or white dwarf 
scenario (although a long observation with XMM-Newton is 
needed to confirm it). In principle, optical spectroscopic ob- 
servations of HD 1 19682 could unveil the presence of a com- 
pact object by means of a radial velo city curve, as has been 
done in the the case of 7 C as itself (iHarmanec et al.ll2000l : 
iMiroshnichenko et al.ll2002h . However, Be stars in X-ray bi- 
naries typically have orbital periods in the range 10-200 d, 
so a long-lasting observing program would be needed. In 
addition, the number and intensity of emission lines along 
the whole spectrum would render difficult to obtain accu- 
rate radial velocity measurements. Certainly, to obtain a ra- 
dial velocity curve is not a straightforward task in the case of 
HD 119682. 

Whatever the nature of HD 1 19682/1 WGA J1346.5-6255, 
we stress that it is not related to the background SNR 
G309.2-00.6. 

7. SUMMARY 

After a study of die sources IWGA J1346. 5-6255 and 
SNR G309.2-00.6 with Chandra and XMM-Newton and of 
HD 119682 with the NTT, we have obtained the following 
results: 

• The IWGA J1346. 5-6255 source is consistent with a 
point source. No jets have been observed with Chan- 
dra. We derive an upper limit on the unabsorbed 
flux from any unseen jets of 2.6x10"'^ ergs cm"^ s"' 
(0.5-7.5 keV), which corresponds to a luminosity of 
5.3 X 10"" ergs s"' at a distance of 1.3 kpc. 

• The X-ray spectrum of IWGA J1346.5-6255 is best 
fitted with two-component models including a thermal 
component. A two-component power-lawH-MEKAL 
model provides the best fits to the Chandra and XMM 
data. The hard component dominates the X-ray emis- 
sion. 

• The Chandra and XMM-derived A^h of (1- 
3)xlO^' cm"^ (conservative range encompassing 
all models used) is lower than that towards the SNR: 
(4-ll)xl02' cm-2 (3cr range). 



than that observed with XMM-Newton. Variability on 
short timescales (a few 100 seconds) has been also de- 
tected at a level of 45-55%, with a possible period of 
^1500 s in the Chandra and XMM data sets, which 
could be the rotational period of a neutron star or a 
whide dwarf. 

• IWGA J1346. 5-6255 is coincident with the optical star 
HD 119682 whose spectrum is that of a B0.5 Ve star. 
It is a blue straggler in the ^45 Myr old open cluster 
NGC 5281. The distance is 1.3±0.3 kpc, significantly 
lower dian flie distance to the SNR G309.2-00.6. 

• Based on the above optical and X-ray spectral and tim- 
ing properties of IWGA J1346. 5-6255, we conclude 
that it is unlikely to be a microquasar associated with 
the SNR G309. 2-00.6, and that it is most Hkely a new 
member of the growing class of 7-Cas objects. 

• Further X-ray observations are needed to study the vari- 
ability of the source, and further explore the bump 
near 6.7 keV. Such a study will confirm the nature of 
IWGA J1346.5-6255 as a 7-Cas like object. 

• The search for a compact object associated with the 
G309.2-00.6 requires a deep Chandra observation of 
the SNR. 
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Right Ascension (J2000) 

Fig. 1. — Chandra image (inset) overlayed on to a radio image of the supernova remnant G309.2-00.6, spanning 21'x21', obtained with MOST at 0.843 GHz, 
with a resolution of 43" (Whiteoak & Green 1996). 
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13H7"^30' 47^^00^ 46^30" 46"'00' 45"^30^ 

Right Ascension (J2000) 

Fig. 2.— The XMM-Newtonl?fi image of IWGA J1346.5-6255 and G309. 2-00.6, spanning 20' x20', shown in logarithmic scale. The overlaid black contour 
corresponds to the brightest emission from G309.2-00.6 detected at 0.843 GHz with MOST. The overlaid circle denotes the region extracted from the XMM data 
for the spectral analysis of the diffuse emission from the SNR (see §4). 
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Fig. 3. — Top left: Chandra image of IWGA J1346.5— 6255 in tlie 0.3-10 keV energy range. The three other images correspond to the soft (top right), medium 
(bottom left) and hard (bottom right) energy ranges. The concentric circles have radii of 2 (barely visible), 4, and 6 arcseconds. The white line indicates the axis 
of symmetry of the G309.2— 00.6 radio lobes, the images have been smoothed (minimum smoothing scale = 1 pixel, minimum signal-to-noise ratio = 2). There 
is a hint of extension in the soft energy band along the symmetry axis of the radio lobes, although the low count rate indicates that it is most likely background 
noise. 
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Fig . 4. — Radial profile of IWGA J1346.5— 6255 (solid bars with 90% confidence en'ors) and of Chandra's PSF (dashed bars) in the 0.3-10.0 keV range (see 
33.1l for details). The data are consistent with a point-like source. 
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Fig. 5. — The Chandra spectrum of the point source IWGA J1346.5— 6255 fitted with an absorbed power-law model. The ratio of data to fitted model is shown 
in the bottom panel. An excess is seen around 1 keV. 




Fig. 7. — The Chandra spectrum of the point source IWGA J1346.5-6255 fitted with an absorbed power-law+MEKAL model. The dot-dashed line represents 
the contribution from the MEKAL soft component, while the dashed line is for the harder power-law component. The ratio of data to fitted model is shown in the 
bottom panel. 
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Fig. 8. — The XMM-Afewto« PN (green), MOSl (black) andMOS2 (red) spectra of of the point source IWGA J1346.5-6255 fitted with an absorbed power-law 
model. The bottom panel displays the ratio of data to fitted model. 
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Fig. 10.— The XMM-Newton PN (green), MOSl (black) and M0S2 (red) spectra of IWGA J1346.5-6255 fitted with a two-component absorbed MEKAL 
model. The lower dashed lines indicate the contribution from the softer component (see text for details). The bottom panel displays the ratio of data to fitted 
model. 
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Fig. 12. — The XMM-Newton PN (green), MOSl (black) and M0S2 (red) spectra of IWGA J1346.5-6255 fitted with a two-component absorbed power- 
law+MEKAL model. The lower dashed lines indicate the contribution from the softer (MEKAL) component. The bottom panel displays the ratio of data to fitted 
model. 
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Fig. 14. — The normalized IWGA J1346.5— 6255 lightcurve obtained with Chandra after background subtraction and binning with 150 s (data points with Icr 
error bars), smoothed with a running window of 450 s (solid line) and the sinusoidal fit to the binned data (dashed line). The time origin corresponds to 2004 
December 26 at 09:22:20 UT (MJD 53365.3905). The data show significant variability and a quasi-periodic signal of ~1500 s. 




Fig. 15. — The normalized IWGA J1346. 5-6255 lightcurves obtained with XMM-Newton PN (top, green), MOSl (middle, black), and MOS2 (bottom, red) 
data after background subtraction and binning with 150 seconds. Error bars are at the Icr level. The time origin corresponds to 2001 August 28 at 02:52:50 UT 
(MJD 52149.1200). The data show significant variability with an amplitude at the 45-55% level with a timescale of a few 100 seconds (see §4.1 for details). 
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Fig. 16. — Classification spectrum of HD 1 19682, taken with EMMl on the NTT on 2003 June 5th. The most prominent emission lines and those absorption 
lines useful for classification purposes are labeled. For readability, none of the He I transitions typical of B-type spectra has been labeled. The C III 4650 A line 
is not labeled either. Surprisingly, it is clearly weaker than the neighboring O II (+N II) 4642 A line, suggesting a moderate N enhancement and C depletion. A 
classification of B0.5 Ve is inferred. 
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TABLE 1 
Photometry of HD 1 19682. 



Ref. V (mag) B-y(mag) 

Moffat & Vogt (1973) 7.98±0.01 0.13±0.01 

Humphrevs (1975) 8.00 0.20 

Drilline (1991) 8.06±0.01 0.22±0.01 

H0g et al. (2000) 7.90 ±0.01 0.08 ±0.02 

Sanneretal. (2001) 7.91±0.04 0.16±0.06 

Weighted average 7.97 ±0.06 0.16 ±0.05 



No te. — The Tycho-2 magnitudes Ma g et alj 
120001) have been converted from VT and BT to 
V and B-V magnitudes following their prescrip- 
tion. We ha ve assumed an error of 0.05 mag 
for the lHumphrevs (.1975 ) photometry to compute 
the weighted average. Errors are at the Icr level. 
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TABLE 2 

Spectral FIT PARAMETERS of IWGA J1346. 5-6255 obtained 

WITH Chandra. 



Model 


Parameters 




Fitted values 


W axJa yjyjVy d 










r 




1.41 (1.19-1.65) 




n nrm ^ V 1 ^ 




1 n 28-1 961 

X.JO v^x,z.o — 








n 463 C611 




flux V 1 n~^^ pra cm~^ 

llUA \ /\ 1-\J tlgO dll 


s-h 


1 n en sq-i 421 

X.XJ — X ."T^y 


wnV^ Q * in f*Vfl 1 


/ViiCxlf)'^'' cm"'''* 




15 en 10-0 211 




kT (keV) 




26.9 (9.7-79.9) 




nnrm y 1 0"^ ^ 

llv/J-J-lJ. \^ /\ XV/ J 




7 02 ("5 92-9 1 61 








U.H-OU v,ui _^ 




nux A lu ergs ^ni 


s ) 




wabs * (mekal+mekal) 


yVHvXllr cm ) 




u.iy (u.iz— u.zo) 








rO 71-1 29^ 

V''* ' — X«i.7_/ 




nornii (xlO~') 




3.10 (1.26-5.66) 








80 (>14) 




norm2 (xlO"*) 




8.2 (5.9-9.4) 








0.342 (59) 




flux (x 10 ergs cm ^ 




1.16(1.10-1.23) 


wabs*(power+mekal) 


iVH(xlO^^ cm"^) 




0.19(0.11-0.32) 




r 




1.24(0.98-1.52) 




normr (x 10^) 




1.30(0.98-1.53) 




kT (keV) 




0.97(0.69-1.34) 




norm^T- (xlO"') 




3.38 (1.09-5.85) 








0.341 (59) 




flux (xlO ergs cm ^ 




1.18 (0.89-1.53) 



Note. — All fits have been performed in the 0.5-7.5 keV en- 
ergy range. The quoted fluxes are unabsorbed in the same energy 
range. The MEKAL models assume solar abundances. AH confi- 
dence ranges are 90%. 
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TABLE 3 

Spectral fit parameters of IWGA J1346. 5-6255 obtained with 

XMM-Newton. 



Model 



Parameters 



Fitted values 



wabs*power 



norm (xlO ) 
flux (xlO ergs cm ^ s 



wabs*mekal Nyi{y.W^ cnT^) 

kT (keV) 
norm (x 10"^) 

flux (x 10"'^ ergs cm"^ s"') 



0.22 (0.20-0.24) 
1.69 (1.64-1.73) 
3.77 (3.54-3.98) 

1.04 (418) 
2.22 (2.20-2.24) 



iVH(xl0^^cm-^) 

kTi (keV) 
normi (xlO"*) 

kT2 (keV) 
norm2 (xlO"') 



0.15(0.14-0.16) 
9.2 (8.5-10.8) 
1.15(1.11-1.19) 

1.12(418) 
2.15(2.0-2.18) 



wabs*(mekal+mekal) 



flux (X 10*12 



ergs cm s 



normr (x 10 
kT (keV) 
normtr (xlO*^) 

xl M 

flux (X 10 '2 ergs cm ^ s ') 



0.16(0.15-0.17) 

1.7 (1.4-2.0) 
1.07 (0.42-3.49) 
13.0 (10.6-15.6) 
1.08(0.9-1.12) 

1.025 (416) 
2.16(2.15-2.17) 



wabs*(power+mekal) 



JVH(xl0^^cm-^) 

r 

norm (X lO""*) 
fiiine (keV) 
(Tii„e (keV) 
normiine (x 10^) 

xl (y) 

: 10 '2 ergs cm ^ ! 



0.21 (0.19-0.24) 
1.60(1.53-1.67) 
3.3 (1.6-3.0) 
1.4(1.1-2.2) 
4.6 (2.5-16.8) 
1.005 (416) 
2.20 (2.18-2.22) 



wabs*(power+Gaussian) 



flux(: 



0.23 (0.21-0.24) 
1.71 (1.66-1.76) 
3.85 (3.63-3.97) 
6.77 (6.65-6.91) 
0.19(6x10-3-0.32) 
4.33 (2.39-6.87) 

1.012(415) 
2.20(1.88-2.26) 



Note. — All fits have been performed in the 0.5-8.5 keV energy 
range using PN, MOSl and MOS2 data. The quoted fluxes are unab- 
sorbed in the same energy range. The MEKAL models assume solar 
abundances. All confidence ranges are 90%. 



